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Abstract 

We consider the baryon octet and decuplet magnetic moments in a simultaneous expansion 
in m s and 1/N taking Np/N ~ 1, where N is the number of QCD colors and Np is the 
number of light quark flavors. At leading order in this expansion, the magnetic moments 
obey the non-relativistic quark-model relations. We compute corrections to these relations 

using an effective lagrangian formalism which respects chiral symmetry to all orders in 

l li 

the 1/N expansion. Including corrections up to order m s , we find 8 relations among 
the 9 measured octet and decuplet magnetic moments; including corrections up to order 
1/N and m s , we find 4 remaining relations. The relations work well, and suggest that 
the expansion is under control. We give predictions for the unmeasured decuplet magnetic 
moments. 
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1. Introduction 



In this paper, we consider the baryon magnetic moments in a simultaneous expansion 
in m s and 1/N, where N is the number of QCD colors. There has been a recent revival of 
interest in the 1/N expansion for baryons, started by the results of ref. [1]. For example, 
it was shown that many interesting large- N relations have corrections starting at ~ 1/N 2 . 
These results have been extended using a number of different methods [2] [3] [4] [5] . We will 
use the formalism of ref. [4] . This formalism is based on an exact relativistic treatment of 
QCD, yet makes direct contact with the static quark model. On a more practical level, it 
allows us to write an explicit effective lagrangian in which chiral symmetry is kept manifest 
to all orders in 1/N. 

When the number of light flavors Np > 2, the SU(Np) flavor representations of 
baryons grow with N, and so there are ambiguities in how to extrapolate physical baryon 
states to N > 3. We use the approach of ref. [5], where it is shown that the 1/N expansion 
can be formulated to include all of the states in the SU(Np) flavor representations for 
arbitrary N and N F . The physical results for N = 3 can be obtained without having to 
identify the physical baryon states with particular states for N > 3, and SU(N F ) flavor 
symmetry is also kept manifest in this approach. This expansion is well-defined even if 
Np/N ~ 1 [5]; we will work in this limit, since Np = N = 3 in the real world. 

In the large-iV limit, the magnetic moments obey the non-relativistic quark model 
relations [6]. We find that the leading corrections to these results are suppressed relative 

1 /2 

to the leading terms by order 1/N, m s , and m 8 . We will assume that 0(m s ) and 
0(1/ N) corrections are both 0(e) ~ 30%, and carry out the expansion consistently to 
0(e). Including the 0(m l J 2 ) = 0(e 1//2 ) corrections, we find 8 relations among the 9 
measured magnetic moments for the octet and decuplet baryons; including corrections up 
to 0(e), we find 4 surviving relations. These relations agree well with data, providing 
evidence that the combined 1/N and chiral expansions work well for baryons. We then 
predict the unmeasured decuplet magnetic moments including all corrections up to 0(e); 
these may be measured at CEBAF. 

2. Formalism 

In this section, we briefly review the formalism we will use to obtain our results. 
In ref. [4], it was shown how to write an effective lagrangian describing the low-energy 
interactions of baryons with the pseudo-Nambu-Goldstone bosons (PNGB's) for large N. 
The PNGB's are described in the standard way: the field 

i(x) = e m ^/f, (1) 
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is taken to transform under SU(N F ) L x SU(N F ) R as 



(2) 



where this equation implicitly defines U as a function of L, i?, and £. For Np = 3, the 
meson fields are 



v/2 



V 



V2 n T vV 
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TV 
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(3) 



Note that the rj' is not light if N F /N ~ 1, and is therefore not included. The effective 
lagrangian is most conveniently written in terms of the hermitian fields 



which transform under SU(Nf)l x SU(Nf)r as 



(4) 



(5) 



We can incorporate SU{Np) breaking due to m s ^ by including the quark mass 
spur ion (for arbitrary Np) 



m. 



m s S, S 



V 



\ 



(6) 



transforming under SU(N F ) L x SU(N F ) R as m g 
define the even-parity field 



Lm q R} . We find it convenient to 



m = — (^m q £ + h.c.) I— > UmU^ . 



In this notation, the leading terms giving rise to meson interactions are 

£ = / 2 tr(A^)+/3/ 3 tr(m) + ---, 



(7) 



(8) 



where (3 is a coupling and / = ~ 93 MeV; / ~ viV in the large- N limit [7]. 

We now discuss the baryon fields. Because the baryon mass is of order ./VAqcd, we 
can describe the baryons using a heavy-particle effective field theory [8]. We write the 
baryon momentum as P = M v + k, where M ~ N is a baryon mass and v is a 4- velocity 
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(v 2 = 1) which defines the baryon rest frame. We then write an effective field theory in 
terms of baryon fields whose momentum modes are the residual momenta k. This effective 
field theory gives an expansion in 1/Mq around the static limit. 

For N large, the baryon SU(Nf) representations are large, and it is convenient to use 
a compact notation to keep track of baryon flavor quantum numbers. We use a Fock-space 
notation in which the baryons fields are written 

\B(x)) = B a ^- a » a » (x) «t iai • • • at NCCN |0). (9) 

The a^'s are bosonic creation operators which create a "quark" with definite flavor and 
spin, and |0) is the Fock "vacuum" state; ai,...,a,N are SU(Np) flavor indices and 
ai,...,ct7v=|,| are spin indices in the rest frame defined by v. 

Under SU(Nf)l x SU(Nf)r, the baryon fields transform as 

^gaiai---ajvajv ^ JJ a l ^ .JJ a N ^ j^b 1 a 1 ---b N ct N (10) 

where U is defined in eq. (2). B " transforms under a highly reducible representation of 
SU(N F ) : but we will have to carry out calculations explicitly only for N = N F = 3. 

Meson-baryon interactions are written in terms of operators constructed from the 
creation and annihilation operators. For example, the leading terms involving baryons can 
be written 

C = (B|w"V M |B) + g (B\{A^}\B) + ■■■, (11) 
where cr M = (^7 M — "1^)75 is the spin matrix and we define 

{A»a,} = aUAT b M a ^, (12) 

etc. The chiral covariant derivative acting on the baryon fields is defined by 

V^\B) = {d^-i{V^})\B). (13) 

The coupling g can be determined from matrix elements of the AS = 1 axial current 
measured in semileptonic hyperon decays. We obtain g = 0.83 ± 0.08, where the error is 
obtained by assigning a 30% uncertainty to the higher-order corrections. This value should 
be used with caution, since the 5 , L r (A'"i?)-breaking corrections are known to be large [8] [9]. 

In this notation, the leading iV dependence of an arbitrary term in the effective la- 
grangian is given by associating a factor of l/iV r-1 with every r-body operator (that is, an 
operator constructed from r creation and r annihilation operators), and a factor of 1/N 
for every explicit flavor trace [4]. The reason for these rules is that an r-body operator 
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can arise only from quark-level diagrams involving at least r — 1 gluon exchanges, and 
flavor traces arise from quark loops. Each gluon exchange or quark loop gives rise to a 
suppression of 1/A, yielding the rules given above. For more details, see ref. [4]. According 
to these rules, the coupling g in eq. (11) is order 1 in the large- A?" limit. 

3. Magnetic Moments 

We now apply the formalism discussed in the previous section to the magnetic mo- 
ments. 

3.1. Leading Order 

At leading order in the large-A and chiral limits, the baryon magnetic moments are 
described by a single term in the effective lagrangian: 

8£ = ^v^(B\{Q*„}\B), (14) 
where W v = ±e^ Xp F Xp , e i 2 3 = +1, and 

Q = \ (£ f QU ] + ZQrO ^uquI (15) 

Here (at N F = 3) 



Ql = Qr = 




are the left- and right-handed quark charge spurions. The parameter A ~ 1 GeV is the 
chiral expansion scale; naive dimensional analysis leads us to expect that ao ~ 1. The 
magnetic moments arising from eq. (14) are given by 

»B>B^B = -^(B>\{Q*i}\B), (17) 

where cr B > B is the matrix element of the spin matrix between the states \B) and \B') 
normalized so that its maximal value is +1. (We use non-script capital letters to refer 
to specific baryon states.) The operator {Qcx 3 } has matrix elements O(A), so that the 
leading contributions to the magnetic moments are O(N).* 

* We could normalize the electric charge of the quarks to be of order 1/A so that the 
baryons have electric charge of order 1 in the large-A limit. If we did this, the magnetic 
moments would be 0(1) in the large-A limit. Such a rescaling would modify the formulas 
which follow in a trivial way, and would not affect our results for A = 3. 
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At this order, there are 9 measured octet and decuplet baryon magnetic moments 
determined by a single unknown constant a / A. There are therefore 8 relations among the 
magnetic moments. Of these, 6 are the Coleman-Glashow relations [10] 





= o, 


(15%) 


+ Hp + 


= o, 


(10%) 


\l n — 2flA 


= o, 


(40%) 




= 0, 


(55%) 




= 0, 


(40%) 




= 0, 


(5%) 



(18) 



which hold in the limit of exact SU(3) independent of the 1/N expansion; the remaining 
relations can be taken to be the quark-model relations 

3/x n + 2fi p = 0, (3%) (19) 
Vta- + l* P = 0, (35%) (20) 

which are consequences of the large- N limit. The numerical accuracy indicated is defined 
by dividing the numerical value by the average of the positive and negative terms on the 
left-hand side. If we perform a best fit, the average deviation is 0.3 /i N . 

3.2. 1/N Corrections 
The 0(1/N) corrections to the magnetic moments arise from the term 

8C = ^v,F^{B\{Q}{a v }\B). (21) 

Because matrix elements of the operator {<2}{o" M } can be O(N), the contributions to the 
magnetic moments arising from this term are 0(1) in the large- N limit. Including this 
term, the Coleman-Glashow relations in eq. (18) still hold, but the quark-model relations 
eqs. (19) and (20) no longer hold. Note that the relation eq. (19) involves states with zero 
strangeness, and therefore receives further corrections only from isospin breaking, which 
are expected to be about 5%. The experimental deviation of this relation is therefore a 
direct measure of the 1/N corrections, and we have no understanding of why the these 
corrections are so small in this case. 
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3.3. SU(Np) -breaking Corrections 



In chiral perturbation theory, the leading SU(Np) -breaking corrections generally arise 
from loop graphs with PNGB intermediate states. Such graphs can give nonanalytic depen- 
dence on the quark masses. In the case of the magnetic moments, the leading dependence 

1 /2 

on m s is ~ m s and ~ m s \nm s [11]. The diagrams which give rise to these corrections 
are shown in fig. 1. 

l li 

The graph which gives rise to the ~ m s corrections is easily evaluated using the 
meson-baryon coupling from eq. (11): 

fig. la = (B'\{T A a»}{T B o- x }\B) tr(T A [T B , Q}) 

1 (22) 
d 4 k (2k + q) fI K(k + q) x 1 ' 



^B'B- J B'B = g -T^( B '\°K\B\ (24) 



(2tt) 4 (k 2 -M\ + iO+)((k + q) 2 - Mf + iO+)(k ■ v + i0+) ' 
where the T A are SU(N F ) generators normalized so that tx(T A T B ) = b~ AB , and 

M 2 AB =^tr({T A ,T B }m q ) (23) 

is the mass-squared matrix of the PNGB's (see eq. (8)). We have neglected the 0(m s ) mass 
differences between baryons in the same SU(Np) multiplet, as well as the 0(1/N) mass 
differences between octet and decuplet baryons. (Including these effects gives corrections 
suppressed by m s and/or 1/N.) The sum over all intermediate spin states is then included 
in eq. (22), with the large- N relations properly taken into account. Evaluating eq. (22) 
gives rise to a contribution to the magnetic moments 

9 2 eM K 

where 

O j K = (N + N F ) tr[Q(l - S)} {Sai} _ {S}{Qa^} + {Q}{S<r'}. (25) 

This gives contributions to the magnetic moments which are 0(Nml /2 ) = OiNe 1 / 2 ). 

Including these contributions along with the leading term in eq. (14), we obtain 8 
relations valid to 0(Ne 1 ^ 2 ) which are independent of g. One of these is the quark-model 
relation eq. (19). The remaining relations can be written 

A*E" + Hn + Hp 

H~.o — 2(i E - 

Ht,+ - 2v / 3 / u E o A + ii p 

1 „ (26) 

HZ+ + ~^HZ°A + HA + Hp 

Hz~ + HY.+ - A»e- - Hp 
Hn- - 1*3° - He- 



o, 


(10%) 


o, 


(4%) 
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(7%) 
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(2%) 


o, 


(4%) 


0. 


(2%) 
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In addition, there is one relation which depends on g, 



^--^-= 9 -^jf ("0.51 = -2.0). (27) 



(We use / = fx — 114 MeV in the evaluation.) The relations which are independent of g 

1 /2 

work much better than the (/-dependent relations: a fit including the 0{Nm s ) corrections 
and treating g as a free parameter has an average deviation of 0.08 n N . The nonanalytic 
corrections have the right sign, but their predicted magnitude for the lowest-order value 
g ~ 0.8 is too large. However, we expect that including the SU(Np>) -breaking corrections 
in the fit to the semi-leptonic decays will substantially decrease g [8] [9], and it is not clear 
to us that the large discrepancy in eq. (27) indicates a breakdown of the expansion. 

The vertex graph (fig. lb) and wavefunction graphs (fig. lc) combine to give the 
contribution 

figs- lb,c = ^e^V [{T A a a }, [{T A ^h{Qa"}]]\B) 



X 



(2tt) 4 (k 2 -M% + i0+)(k • v + i0+)((k - q) ■ v + i0+) 



The double commutator can be written as a 1-body operator; this contribution is therefore 
~ 1/N times a 1-body operator, and can be at most 0(m s lnm s ) = 0(Ne 2 \ne). This is 
negligible compared to the O(Ne) counterterms which we will discuss below. The graph in 
fig. Id gives a contribution which is 1/N times a 1-body operator, which is negligible for 
the same reason as for eq. (28); the graph in fig. le gives no contribution to the magnetic 
moments 



The leading S 1 ^ (iVi?)-violating counterterms are 



5£ = ev^ v {B\ 



|B). (29) 



These counterterms give 0(Nm s ) = 0{Ne) contributions to the magnetic moments. There 
are also 0(Nm s ) contributions from the loop graph in fig. la, but these have the same spin- 
flavor dependence as the counterterms considered above, so we will not need to evaluate 
them explicitly. Including the counterterms in eq. (29) and the nonanalytic corrections 
computed above, we obtain the relations 

^=0 + 2^ E + + 2^ E - + fi n = 0, (10%) 
(is- + 4// E - + 2v / 3^ e oa + + 8^ n = 0, (2%) 

you ) 

Ha- + 4^ H o - 3/x H - + 8^i S + + 5^- - 3fi p + [i n = 0, (8%) 
4^ s o - - |it E - + 4v / 3 / u E o A - 6fi A + 4/x n = 0, (6%). 
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The last of these relations was noted in ref. [12], and is valid to 0(m s ) independent of the 
1/N expansion. The average deviation of a fit which treats all of the counterterm couplings 

1 /2 

as free parameters is 0.08 fi Nl the same as the fit including only the 0(Nm s ) terms. 



4. Predictions 

We can use the results obtained above to predict the unmeasured decuplet magnetic 
moments including all contributions up to O(Ne). These may be measured in the future 
at CEBAF. 

Up until now, the fits were used only to give a rough idea of how well the expansion 
works. We now use the fit to the magnetic moment data to make predictions, so we give 
some details on the fit and the treatment of errors. We add a theoretical uncertainty of 
0.1 in quadrature with the experimental error to obtain the error on the individual 
magnetic moments, and then perform a x 2 fit. This theoretical uncertainty is approxi- 
mately the average deviation of the fit when corrections up to O(Ne) are included, and 
is consistent with the expectations of dimensional analysis: the largest contributions not 
included are mj 2 nonanalytic terms 

M 3 M 

6 ^Y^p^^-0.2^ N . (31) 
We then obtain the predictions (in Bohr magnetons) 





= 5.9 ±0.4, 




= 2.9 ±0.2, 




= 3.3 ±0.2, 




= 0.3 ±0.1, 


A*E*~ 


= -2.8 ±0.3, 




= 0.65 ±0.2, 




= -2.3 ±0.15, 



(32) 



The error quoted is purely the formal error from the fit. Note that in the limit of exact 
SU (3) flavor symmetry, the decuplet magnetic moments are proportional to their charges. 
Eq. (32) gives definite predictions for the pattern of SU(3) violation which would be very 
interesting to test. 

We make no prediction for the A magnetic moment. The leading contribution to 
this magnetic moment comes from terms such as 

6C= irh tr(mS) ^Ap^ A (siK}|£), (33) 
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where 6 4 ~ 1. (The reason for the factor of 1/N in the coefficient is that the quark 
diagrams which contribute to this term contain a fermion loop; see refs. [4] [5].) This gives 
a contribution to the magnetic moments 0(1/N) = 0(Ne 2 ), which is higher order than 
the terms we are keeping. 

Decuplet magnetic moments were also considered in a chiral expansion in ref. [13]. 

1 /2 

This paper computes 0{m s ) and 0(m s lnm s ) 5 , L r (3)-violating contributions to the de- 
cuplet magnetic moments without making use of the 1/N expansion. Their predictions 
differ considerably from ours: for example, their predicted values for and fx-.* are 
significantly smaller than ours. It is worth noting that they do not include 0(m s ) countert- 
erm contributions, which are not expected to be significantly smaller than the 0(m s lnm s ) 
which they compute, whereas the largest terms which are omitted in our analysis are sup- 
pressed by a power of the expansion parameter. We hope that experiment will be able to 
decide between these different predictions in the future. 

The radiative decays of decuplet baryons also receive contributions from the same 
terms that give rise to the magnetic moments. However, the momentum transfer for the 
process T — > S7 is 

|g1 ~ M T - M B (34) 

and the loop graph in fig. la has nontrivial dependence on \q\/ (Mt — Mb) ~ 1. This makes 
the coefficients of the effective operators in the expansion different from the ones appearing 
in the expansion of the magnetic moments. In principle, this difference is computable, but 
it depends on the value of the axial coupling g, which we know is not well determined. We 
will return to these issues in a future publication. 
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Figure Captions 

Fig. 1. Feynman graphs giving rise to nonanalytic corrections to the baryon magnetic 
moments. Fig. lc is wavefunction renormalization; fig. le does not contribute to the 
magnetic moments. 
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